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Abstract 

In this paper we prove an optimal error estimate for the H(curl)-conforming projection 
based p-interpolation operator introduced in [L. Demkowicz and I. Babuska, p interpolation 
error estimates for edge finite elements of variable order in two dimensions, SIAM J. Numer. 
AnaL, 41 (2003), pp. 1195-1208]. This result is proved on the reference element (either 
triangle or square) K for regular vector fields in H''(curl, iiT) with arbitrary r > 0. The 
formulation of the result in the H(div)-conforming setting, which is relevant for the analysis 
of high-order boundary element approximations for Maxwell's equations, is provided as well. 
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1 Introduction 

This paper concerns the H(curl)-conforming interpolation of regular vector fields by high or- 
der polynomials on the reference triangle or square and the corresponding interpolation error 
estimation. To the best of our knowledge, the first paper related to this subject is the paper 
by Suri from 1990 [23j. In that paper error estimates (in terms of both the mesh parameter 
h and the polynomial degree p) were derived for classical Raviart-Thomas (RT) and Brezzi- 
Douglas-Marini (BDM) interpolation operators on the reference square Q (note that in 2D 
these H(div)-conforming so-called face elements and the H(curl)-conforming edge elements of 
the Nedelec type are isomorphic). The estimates obtained were not optimal with respect to p 
and later, in [22], they were improved to e-suboptimal p-estimates for sufficiently regular vec- 
tor fields u (namely, for u G H'^'(div,Q) with r > ^). These results were further extended by 
Ainsworth and Pinchedez in [Ij (to meshes with hanging nodes, weighted Sobolev regularity 
of approximated functions, exponential convergence on graded meshes) and in [2] (to Brezzi- 
Douglas-Fortin-Marini elements with non- uniform distribution of polynomial degrees). In 3D 
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the corresponding results were obtained by Monk in [19] and by Ben Belgacem and Bernardi in 
[7]. All mentioned papers deal with quadrilateral or hexahedral elements and the proofs therein 
essentially rely on expansions in terms of orthogonal (namely, Legendre) polynomials. An ap- 
plication of this approach to triangular or tetrahedral elements does not seem to be feasible. 
Another drawback of classical interpolation operators for edge (or face) elements is the lack of 
stability (with respect to p) for low-regular fields. 

A breakthrough in H(curl)-conforming p-interpolation analysis, i.e., the construction of an 
interpolation operator which works equally well on both triangular and quadrilateral elements 
and also for low-regular fields, was achieved relatively recently by Demkowicz and Babuska 
in [T3]. These authors have introduced and analyzed H^- and H(curl)-conforming projection- 
based p- interpolation operators satisfying the commuting diagram property (de Rham diagram). 
This property and the corresponding p-interpolation error estimates have immediate applications 
to the analysis of high-order finite element (FE) discretizations of the time-harmonic Maxwell's 
equations. In particular, they are critical to prove the discrete compactness property (which 
implies the convergence of FE approximations for Maxwell's equations) and also useful for the 
error analysis (see jlH llOj). Moreover, the interpolation operators in were constructed 
to allow polynomial degrees to vary from one element to another. This has been done by 
assigning to each element an "internal" polynomial degree and a sequence of (possibly lower) 
"edge" degrees. Such a construction of interpolation operators is essential for the analysis of 
exponentially convergent /ip-approximations and /ip-adaptive schemes. In [jl5j these results have 
been extended to the 3D case. 

The error estimates presented in |14j for both p-interpolation operators are suboptimal. 
Furthermore, the error estimate of the H(curl)-conforming interpolation is available only for 
low-regular vector fields in H^(curl, K) with r E (0, 1) (here K is either the reference triangle or 
square). Though these drawbacks are not essential for the convergence analysis of FE approxi- 
mations, the corresponding improvements would be advantageous for the error analysis. In this 
paper we show that an optimal estimate for the error of the i/^-conforming p-interpolation can 
be obtained in the if^-semi-norm. Using this result we then prove an optimal error estimate for 
the H(curl)-conforming p-interpolation operator applied to a vector field u G H'"(curl, JT) with 
arbitrary r > 0. In the proof we rely on a regular splitting of u € il'''{cml,K) into a curl- free 
component and a complementary vector field of extra smoothness. This splitting is possible to 
recent results related to the regularized Poincare-type integral operators in [13J. 

The paper is organized as follows. In the next section we introduce necessary notation and 
formulate some auxiliary results. In Section [3] we quote [13] to briefly sketch the definition and 
properties of the H^- and H(curl)-conforming projection-based ^-interpolation operators. Opti- 
mal error estimates for both interpolation operators are proved in Section U] (see Theorems 14.11 
and 14. 2p . The paper is concluded with Section [5] where we mention some simple extensions of 
our results including the /f(div)-conforming p-interpolation and /ip-estimates. 

Throughout the paper, C denotes a generic positive constant which is independent of p and 
involved functions. 
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2 Notation and auxiliary results 



We will present all technical details only for the equilateral reference triangle T = {x2 > 0, X2 < 
V2,{xi + 1), X2 < (1 - xi)^/3}. The case of the reference square 

Q = (-1, 1)2 (for which the 

arguments are essentially the same) is briefly discussed in Section [5j A generic side of the triangle 
T will be denoted by 

We will use the standard definitions for the Sobolev spaces H"^ (r > 0) of scalar functions 
on the interval I = (—1, 1) and on the triangle T (see, e.g., |18j). The norms in these spaces are 
denoted by || • ||_H-r(/) and || • ||h' (t)) respectively. On the interval / we will also need the Sobolev 
spaces H'^(I) for r G (0,1) which are defined by interpolation. We use the real K-method of 
interpolation (see [18j) to define 

H^{I) = H'.il)) ^ (1/2 < t < 1, < r < t). 

Here, H^q{I) (0 < i < 1) is the completion of C(f>(/) in H\l) and we identify H^{I) and H^{I). 
Note that the Sobolev spaces also satisfy the interpolation property, e.g., 

H^{I) = [l\I),H\I))^^^ (0<r<l) 

with equivalent norms. Furthermore, the Sobolev spaces H^' and on any edge i C dT are 
defined by using the definitions of the corresponding spaces on the interval /. 

Throughout the paper, we use boldface symbols for vector fields. The spaces (or sets) of 
vector fields are denoted in boldface as well (e.g., H''(r) = {H''{T)f), with their norms and inner 
products being defined component-wise. The standard notation will be used also for differential 
operators V = (d/dxi, 8/8x2), div = V •, curl = Vx, and for the Laplace operator A = div V. 

The L^-inner product and the corresponding L^-norm on T are denoted by (•, •)o,t and || • ||o,T) 
respectively. We will use the semi-norm in H^{T) which is defined as 

\u\h^{T) = ||Vn||o,T- 

Furthermore, we will use the space 

H'^(curl,r) := {u G H'^(r); curlu G H''{T)], r > 0, 

and its analog H''(div,r) in the div-setting. In both cases the spaces are equipped with their 
graph norms. For r = we drop the superscript in the above notations: H''(curl, T) = H(curl, T) 
and HO(div,r) = H(div,T). 

We will also need the space H^/'^{8T) which can be defined as the trace space of H^{T) on 
8T with norm 

\W\\m/^(dT) = ,,inf \\U\\m{T). 

Let us introduce the needed polynomial sets. By Vp{I) we denote the set of polynomials of 
degree < pon the interval /, and Vpi^I) denotes the subset of T'p{I) which consists of polynomials 
vanishing at the end points of /. Li particular, these two sets will be used for the edges £ C 8T. 
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Further, Vp{T) is the set of polynomials on T of total degree < p. The corresponding 
set of polynomial (scalar) bubble functions on T is denoted by Vp{T). When considering the 
reference square Q we will denote by Vp{Q) the set of polynomials of degree < p in each variable 
separately. Finally, Vp^'^{T) denotes the polynomial set associated with the second Nedelec 
family of edge elements on T, i.e., Vp'''^{T) = {Vp{T)f (see [H]). The subset of ^^^^(r) 
which consists of vector- valued polynomials with vanishing tangential trace on the boundary dT 
(vector bubble- functions) will be denoted by 'P^^'^'^{T). 

To simplify the presentation we have assigned only one ( "internal" ) polynomial degree to the 
reference element. Following [14] the results extend to the polynomial spaces on T with separate 
polynomial degrees (of possibly lower order) assigned also to traces on the edges of dT. 

When proving interpolation error estimates in Section 0] we will need some auxiliary results, 
which are collected in the next three sub-sections. 

2.1 Polynomial extensions from the boundary 

The problem of polynomial extension from the boundary dT can be formulated as follows: given 
a continuous function / defined on dT such that f\i G Vp{t) for any (. C (9T, find a polynomial 
<I> G 'Pp{T) satisfying ^\qt = f ■ The existence of such an extension, which is stable (with respect 
to p) as a mapping H^/^{dT) H\T), has been proved in [H Theorem 7.4]. In general (i.e., 
for p > 2), the extension is not uniquely defined. To ensure the uniqueness, we will search for 
discrete harmonic extensions: given a continuous piecewise polynomial / of degree p on each 
side i C dT, find a polynomial F := £pf £ Vp{T) such that F = f on dT and 

(VF,V(^)o,T = V(^GP°(T). (2.1) 
Then for any polynomial extension $ G Vp{T) such that ^\dT = / there holds 

I*^Ihi(T) = K*^ ~ ^pf) + '^p/l/i'i(T) = 1*^* — '^p/l/i'i(T) + I'^p/Ihi(T) — l^p/l/i'i(T)- 

Hence, applying the mentioned result of [3], we have 

\£pf\m(T) < C' 11/ 11^1/2 (aT) (2-2) 
with a positive constant C independent of p and /. 

Remark 2.1 (i) The same result as above holds for the reference square Q (see ^ Theo- 
rem 7.5]). 

(ii) If f is a continuous piecewise linear function on dT and F := Sif £ Vi{T), then F = Spf 
for any p > 1. Indeed, F £ Vi{T) C Vp{T), F = f on dT, and for every ip £ Vp{T) there holds 

iVF,Vip)o,T = -(AF,(^)o,T = 0. 
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2.2 Polynomial approximation of scalar functions 

In this sub-section we provide two ^-approximation results for scalar functions with Sobolev 
regularity in one and two dimensions. The first result is a careful revision of [6l Theorem 3.2], 
where an additional log^^^p-term in the error bound appears. 

Lemma 2.1 Let I = (—1,1) and f G H''{I) with some r > ^. Then there exists a polynomial 
fp G Vpil) such that /p(±l) = /(±1) and 

11/ - < Cp-^'^-i/^) ||/||^.(,). (2.3) 

Proof. We outline a proof following the ideas of O |6]. Let / G H^{I), r > 1/2. Then / can 
be expanded as 

oo 

f{x) = ^airj(x), 

i=0 

where Ti{x) is the Chebyshev polynomial of degree i on /. 
Setting 

p 

Ppf{x) ■.= ^aiT,{x)£Vp{I), 
i=0 

one has for any s G (^,min{l,r}) (see [6l Theorem 3.2, Remark 3.2]) 

\\f-Ppf\\HsiI)<Cp-(^-^U\f\\Hr^jy (2.4) 
Moreover, it is easy to show that (j2.4p holds for s = 0: 

oo 

kip 

i=p+l 
-2r\ 



PpfWh^j) < Il(/-^'p/)(cose)||i2(o,^)=^ |a 



< Cp-^'- 5^ |a,|2(l + i2)'-<C7p- uJUHr(iy 

i=p+l 



Analogously, applying the Schwarz inequality we prove that 

2 



\{f - Ppf){±l)\' < ( ^ \aA 

^ i=p+i ' 



oo 



i^Y 



< (2.5) 
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Now we will adjust the polynomial Ppf at the end points of /. First, using the same idea as in 
the two-dimensional case (see [SJ pp. 759-760]), we find two polynomials , V'p £ "Ppil) such 
that 

^^-(-1) = !, ^p-(l) = 0, ^+(-l) = 0, V'+(l) = l, 

and 

U^\\hs^i)<Cp'-'/^ se{0,l). (2.6) 

Then we set 

fp{x) := Ppfix) + (/ - Pp/)(-l) Vp-(x) + (/ - Pp/)(1) V;(x), xel. 

It is easy to check that fp £ Vp{I) and /p(ibl) = /(itl). Furthermore, making use of ()2.4p - ()2.6p . 
we obtain for any s G {0} U (^,min{l,r}) 

\\f - fp\\H-{I) - \\f - fp\\H-{I) 

< 11/ - PpfWH^H) + max |(/ - Ppf){x)\ [UpWhsU) + U^WrsH)) 

< Cp-^''-'^f\\Hrii). 

The inequality in (|2.3p now follows via interpolation between H^{I) and H^{I) for some s £ 
(i,min{l,r}). □ 

In 2D the following approximation result holds (see [4^ Lemma 4.1]). 

Lemma 2.2 Let K he the reference triangle or square. Then there exists a family of operators 
{vTp}, p = 1, 2, . . . , TTp : H^{K) Vp{K) such that for any f G H^{K), k > there holds 

11/ - VIIh^w < Cp-^'-'^WfWH^^K), o<s<k. 

Moreover, iTp preserves polynomials of degree p, i.e., -Kpf = f if f £ Vp{K). 
2.3 The regularized Poincare integral operators 

In [13], Costabel and Mcintosh studied a regularized version of the Poincare-type integral opera- 
tor acting on differential forms in M". They proved, in particular, that this operator is bounded 
on a wide range of functional spaces including the whole scale of Sobolev spaces H^(fl) (r G M) 
on a bounded Lipschitz domain J7 which is starlike with respect to an open ball. Moreover, 
the essential property of the classical Poincare map to preserve polynomials is retained by its 
regularized version. Thus, the results of [13] have immediate applications to the analysis of 
high-order edge elements (see, e.g., [13 [8] and the proof of Theorem 14.21 below) . 

Let us formulate some results of [13] in two particular cases. Namely, we will define two 
Poincare-type integral operators: one operator acts on scalar functions, and the other one acts on 
curl- free vector fields. In both cases the functions and vector fields are defined on the reference 
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element (either triangle or square) K. Denoting by B an open ball in K, let us consider a 
smoothing function 

0eC7°°(M2), supples, je{a)da=l, a =(01,02). 

B 

Then the first regularized Poincare-type integral operator R : C°°{K) — > {C°°{K))'^ (i-e., the 
operator acting on scalar functions) is defined as Rip = {Ri,R2), where 

1 

i?i(x) := — J 6'(a) (x2 — 02) J tip{a. + t{x — a)) dt da, 

B 
1 

i?2(x) := j e{a){xi - ai) j tip{a + t{yi- a))dtda. 

B 

The second operator acting on vector fields is defined as follows: 

2 1 

Au(x) := j 0(a) ^(xj — Oj) j Ui{a + t{yi — a)) dt da, 

B 

where u = (ni,n2). 

The following properties of the operators R and A are easy to check directly (see also [13\ 
Proposition 4.2]): 

(Rl) i? is a right inverse of the curl operator, i.e., 

curl(iiV) = V' yipeH'^iK), r>0; 

(Al) if u is curl-free, then A is a right inverse of the gradient, i.e., 

V(^u)=u Vu G H'' (curio, K) = {u G H''(i^:) ; curlu = in iC}, r > 0. 

Furthermore, the operators R and A satisfy the following continuity properties (see [13\ 
Corollary 3.4]): 

(R2) the mapping R defines a bounded operator H^{K) — > H^+^(i^) for any r > 0; 
(A2) the mapping A defines a bounded operator H^(i^) H^^^(K) for any r > 0. 

We will use the operators R and A to prove the following auxiliary lemma. 



Lemma 2.3 Let u G H'' (curl, ii') , r > 0. T/ien i/iere exist a function G H^^^{K) and a 
vector field v G i¥'^^{K) such that 

u = VV^ + v. (2.7) 

Moreover, 

l|v||H-+i{i^) < C'l|curlu||^r(;i-) and MHr+HK) < C \\u\\iir^Ky (2.8) 

Proof. Since curlu G H'^{K), we use the operator i? to define the vector field v := i?(curlu) G 
W'+'^{K) (see property (R2)). Then 

u = (u - i?(curlu)) + i?(curlu) = (u - i?(curlu)) + v. (2.9) 

The vector field (u — i?(curlu)) G H'^{K) is curl-free due to property (Rl). Therefore, applying 
the operator A to this vector field and using properties (Al) and (A2), we find a function 
tp := A(u — i?(curlu)) G H'^^^{K) such that V tp = u — i?(curlu). Therefore, the decomposition 
of u in (j2.9|) can be written in the form given by (j2.7|) . The inequalities in (|2.8|) are obtained by 
using the continuity properties of the operators R and A: 

\W\\Hr+-^iK) = ll-R(curlu)||Hr+i(i^) < C||curlu||j^r(;^) 

and 

Mh-+^k) = P(u - ^(curlu))||j^r+i(;^) < C (||u||Hr(ii-) + ||i?(curlu)||Hr(A-)) < C \\u\\iir(K)- 

The last inequality relies on specific mapping properties of the scalar curl operator (see, e.g., [S]). 
This finishes the proof. □ 



3 Interpolation operators 

In [T3] two projection-based interpolation operators have been introduced and analyzed. These 
are the i/^-conforming interpolation operator Up : H^^^{T) 'Pp{T) and the H(curl)-conform- 
ing interpolation operator 11™''' : H'"(T) H H(curl, T) 'Pp'^'^(T) (here, r > in both cases). 
Let us briefly sketch the deflnitions of both operators and summarize their properties (see [14j 
for details). 

Let g G H^'^'^{T), r > 0. To define the interpolant Upg, one starts with the standard linear 
interpolation of g at the vertices of T: 

gi G Vi{T), gi = g at each vertex of T. 

Then, for each edge £ C dT, we define a polynomial g2/ by using the projection 

92,1 ^ Vp{t) : 11(5 - 51)1 52/1 1 j^i/2(^) ^ min. (3.1) 
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Extending g2/ by zero onto the remaining part of dT (and keeping its notation), using the 
polynomial extension £.p from the boundary (see Section I2.ip , and summing up over all edges 
we define 

al ■■= E ^p(92,i) e VpiT). (3.2) 

ecdT 

Finally, we define the polynomial bubble by projection in the if^-semi-norm 

g^eV^iT): \{g - gi- g^) - glluif^r) ^min. (3.3) 
Then the interpolant Hp g is defined as the sum 

Ulg:=gi+gP^+gP^erp{T). (3.4) 
Note that, due to the finite dimensionality of gi, there holds for r > 

\\9i\\m+r(T) - X] 1^(^)1 - ^yWcif) < C\\g\\Hi+r(T)- (3-5) 

v. vertices of T 

Now we proceed to the H(curl)-conforming interpolation operator. Given a vector field 
u G W{T) n H(curl,r) with r > 0, the interpolant = H™"^' u G Vp'^'^iT) is also defined as 
the sum of three terms: 

= Ui + U2 + U3. 

Here, ui is the Witney (lowest order) interpolant 



ui = ^ (^j nxudo^ cf)^, 



ecdT 



I 



where n = (721,77,2) denotes the outward normal unit vector to 9T, n x u = 77,17x2 — n^ux with 
u = (7x1,772), and are the standard basis functions (associated with edges €) for V^'^^iT). 
For any edge I C dT one has 



/ 



n X (u — ui) da = 0. 



Hence, there exists a scalar function ip, defined on the boundary dT, such that 

— — = n X (u — ui), 7/; = at all vertices. 
da 

Then, for each edge £, the restriction tpli is projected in the ^-'^/^(£)-norm onto the set of 
polynomials Vp_^_i{i) 

i^2,e e 'Pp+ii^) ■ Me - V'2,£|Ihi/2(q ^ min. 
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Extending 'tp2,e by zero from i onto dT (and keeping its notation) and using the polynomial 
extension £p+i from the boundary we define "02,^+1 '•— ^p+i(V'2/) G "Pp+iiT). Then we set 

"2= E ^.e^p^'m, where 4, = v4,+i. 

ecdT 

The interior interpolant Ug is a vector bubble function that solves the constrained minimization 
problem 

Uf G 7'Ned,0(^) . 

||curl(u - (ui + + uf))||o,T min, 
(u - (ui + + uf), V 0)o,T = V</. G T'0+i(r). 
These interpolation operators satisfy the following properties. 

1°. For r > the operators : U^+''{T) H\T) and n^'^''^ : H''(r) n H(curl,r) ^ 
H(curl, T) are well defined and bounded, with corresponding operator norms independent 
of the polynomial degree p (cf. [14, Propositions 1, 2]). 

2°. The operators 11^ and 11^"'^^ preserve scalar polynomials in Vp{T) and polynomial vector 
fields in V^'"^{T), respectively. 

3°. For r > the following diagram commutes (see Proposition 3 in [14J): 
H^+r{T) ^ H'-(r)nH(curi,r) ^ L2(r) 



n^+1 



TTCUrl 

p 



K-1 (3-6) 



where Hp : Lp'{T) — > VpiT) denotes the standard L^-projection onto the set of polynomials 

4 Interpolation error estimates 

First, we consider the i/^ -conforming p-interpolation operator. We prove that an optimal esti- 
mate can be obtained for the interpolation error measured in the i?^-semi-norm. 

Theorem 4.1 Let g G H^^^'iT), r > 0. Then there exists a positive constant C independent of 
p and g such that 

Iff - ^la\m{T) < Cp^'' Mm+'-iT)- (4-1) 

Proof. Let g G H^^'^{T) with r > 0. Using the operator vTp (see Lemma [2. 2 p and the polynomial 
extension 8p from the boundary, we define the following polynomial bubble function on T: 
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Hereafter, 7tr denotes the standard trace operator, 7tr/ = f\dT- Then, making use of the 
definition of the interpolant Upg (see ()3.2p - ()3.4p ). we have 

|5 - np5l//i(T) < \{9 - 91 - 92) - 93\h^T) <\(.9 - 91 - 92) - ^p\m{T) 
< \9 - ■^p9\m{T) + 91 + 92 - ^p(.ltv{-n-p9)) 



< \\9 - ■^p9\\h^t) + 



m(T) 

£p{ltr{9l) + ^ 92,1 - ltr{T^p9) 

ecdT 



(4.2) 



Here we also used the fact that gi = <?i (7^(51)) = £pi'ytr{9i)) for any p > 1 (see Remark l2.ll 
(ii)). The extension operator £p satisfies (12. 2p . Therefore, applying the continuity property of 
the trace operator jtr ■ H^{T) — > H^^'^{dT), we find 



£pntr{9i)+ X] 52/ - 7tr(vrpfif) 



ecdT 



< \ntr{9l)+ X] 32,£ - Itr {■n-p9) 



ecdT 



< C I 7^(51) + X] ^2,£ - 7^(5) 



ecdT 



- ^ ( II ^ - 7^(5 - 91) 

^ ecdT 



H^/^(dT) 
m/^(dT) 



ltr{9 - T^p9) 

+ ll5-7rp5lbi(T)^- 



//i/2(aT) 



Using this estimate in (14. 2 p we obtain 



Iff - ^,9\miT) <C\\g- -np9\\m{T) + || J] ^2,1 - ltr{9 - 9i) jri/2(ST) 

ecdT ^ ' 

The first term on the right-hand side of (j4.3p is estimated by applying Lemma [ 

\\9 - '^p9\\m{T) < Cp"*" \\g\\H^+r(^T)- 



(4.3) 



(4.4) 



To estimate the second term on the right-hand side of (j4.3p we use localization to the edges 
of the triangle, the definition of the edge interpolants §2/ (see (|3.ip ). and the p-approximation 
result in H^/'^{i) on each edge (see Lemma l2.ip . One has 



X] 92,e - 7tv{9 - 91) 
ecdT 



Hi/2(aT) 



- \\92,e - {9 - 9i)\e\\Hi/2 



ecdT 



ecdT 



9i)\e\\H^/^+'-{e)- 
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Hence, applying the trace theorem for individual edges (see \16\ Theorem 1.4.2]) and estimating 
the norm of gi as in (j3.5p . we prove 

II Yl ^2,^ ~ ^t'-^^ ~ m/vaT\ - ~ 5i||//i+'-(T) < Cp^' \\9\\H^+r(T). (4.5) 

Now the desired error bound in (j4.ip immediately follows from (j4.3p - (j4.5p . □ 

In the following theorem we prove an optimal error estimate for the H(curl)-conforming 
p-interpolation operator 11^"'''. 

Theorem 4.2 Let u E H''(curl, T), r > 0. Then there exists a positive constant C independent 
of p and u such that 

||U - n^"''' u||H{curl,T) < Cp^'' ||u||H'-(curl,T)- (4.6) 

Proof. Given u G H''(curl, T) (r > 0), we use Lemma [2]3] to decompose u as follows: 

u = VV + v, VeF'^+^T), vGH''+i(T). (4.7) 

Moreover (see (jTH]) ). 

I|v||h^+i{T) < ||curlu||j:^r(r) ||V'll_H''-+i(T) < C* ||u||h'-{t)- (4.8) 

Then, applying the interpolation operator 11^'^'^^ and using its commutativity with ^^j^i (see 
(|3.6p ). we write 

n;'^'-' u = n7'(v ^) + n^' v = v(ni+iV) + v. (4.9) 

Since IIp™^ is a bounded operator preserving polynomials (see properties 1°, 2° of 11^™'), one 
has for some fixed e G (0, 1) and for any polynomial Vp G (VpiT))'^: 

||v-n^"'''v||H{curl,T) = l|v- Vp-n™'''(v- Vp)||H{curl,T) 



- ^ ^^'^f, l|v- Vp||h^(7.) + ||curl(v-Vp)||i2(r) 

vpe{Pp(T))^ ^ 

< C inf l|v — VrillHirrv 

vpe{Pp(T))2" ^" ^ ' 

Hence, applying Lemma 12.21 componentwise and using the first inequality in (j4.8p . we estimate 

II V - n^"'' v||H(curl,T) < Cp"" l|v||Hi+'-(r) < C P~'' l|curl u||^r(-jn) . (4.10) 

On the other hand, applying Theorem 14.11 and the second inequality in ()4.8p we obtain 

||V V - V(nl+iV)llH(curl,T) = IV' - Ill+i^\m(T) < Cp-' ||V||hi+'-(T) < Cp-' ||u||h.(T)- (4.11) 

Combining (j4.10p and (|4.1ip we prove (|4.6p by making use of decompositions (j4.7p . ()4.9p and 
the triangle inequality. □ 
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5 Concluding remarks 



The main result of the paper - an optimal error estimate for the H(curl)-conforming projection 
based p-interpolation operator 11™'^^ introduced by Demkowicz and Babuska in [H] - is proved 
for the second Nedelec family of edge elements on the reference triangle T (see Theorem 14. 2p . 
Below we discuss some simple extensions of this result. 

The first Nedelec family of edge elements and the case of the reference square. Using ap- 
propriate polynomial spaces and the same constructions as in Section [31 one can define the 
H(curl)-conforming p- interpolation operator for the first Nedelec family of edge elements in- 
troduced in [2D], and also for both the first and the second families on the reference square 
Q = (—1,1)^. Moreover, properties 1°- 3° of the operator 11^"''^ formulated in Section [3] will 
remain valid in all these cases (see |14[ Section 6]). Thus, the proof of Theorem 14.21 carries over 
without modifications to all the cases mentioned here. 

H (div) -conforming p-interpolation operator. As mentioned in the introduction, the H(curl)- 
conforming p-interpolation operator is critical for the convergence and error analysis of the 
high-order FEM with edge elements for Maxwell's equations in two dimensions. However, when 
a boundary integral formulation of Maxwell's equations (e.g., the electric field integral equation) 
is discretized by the high-order boundary element methods, similar interpolation operators are 
required in the i?(div)-conforming setting for RT- or BDM- surface elements (see [9l|8]). For 
the sake of completeness we will formulate here the main results related to the operator 11™''^ 
in the H(div)-setting. This can be easily done by rotation due to the isomorphism of the curl 
and the div operators in 2D (and, as a consequence, the isomorphism of the first (resp., second) 
Nedelec family of edge elements and the RT- (resp., BDM-) elements). 

Let K be either the equilateral reference triangle T or the reference square Q. We will focus 
the presentation on RT-elements only (in the case of BDM-elements all formulations below are 
essentially the same). Let V^^{K) be the RT -space of order p > 1 on the reference element K 
(see, e.g., [12]), i.e., 

(K) = iv,^iiK)f e ( 2 

Then the H(div)-conforming ^-interpolation operator 11^'^ : H'"(i^) n H(div, /C) — > 7^^'^{K) 
(r > 0) can be defined as the "rotated" H(curl)-conforming interpolation in the same way as in 
Section [3l It satisfies the following properties: 

1°. For r > the operator Uf^ : W{K) n H(div,E:) H(div,E:) is bounded with corre- 
sponding operator norm being independent of the polynomial degree p. 

2°. The operator 11^'^ preserves polynomial vector fields in T^^'^{K). 

3°. For r > the following diagram commutes (cf. (j3.6p ): 



H^+r^K) ^ H''(K)nH(div,i^) ^ L'^{K) 
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where curl = {d/dx2, —d/dxi). 



The following theorem provides an optimal estimate for the interpolation error. 

Theorem 5.1 Let u G H'"(div, r > 0. Then there exists a positive constant C independent 
of p and u such that 

||u - np''u||H(div,/f) < Cp^' l|u||H'-(div,/f)- 

hp-estimates. The projection-based interpolation operators 11^'^''^ and 11^'^ preserve piece- 
wise polynomial vector fields and provide conforming approximations in H(curl) and H(div), 
respectively. Therefore, using the standard Bramble-Hilbert argument and scaling, Theorems l4.2l 
and 15.11 extend to the corresponding optimal /ip-estimates on sequences of quasi-uniform meshes 
of triangles and/or parallelograms satisfying the standard shape regularity assumptions. 
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